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I. INTRODUCTION 



Current knowledge of the quantum structure of space-time suggests a picture very dif- 
ferent from the smooth classical one. Space and time are expected to be fundamentally 
discrete such that both change only in steps. Still, a transition between both pictures must 
be possible in order to understand the emergence of a classical world on large scales from 
the fundamental quantum world. Also for practical purposes such a transition is helpful in 
a semiclassical approximation. A technically and conceptually important question is where 
the classical picture starts to make sense or, when going to smaller and smaller scales, where 
it breaks down. The change of scale can happen either computationally, i.e. by looking at 
smaller and smaller scales in a coarse-grained approximation which can then be used, for 
instance, to understand the breaking or deformation of classical symmetries, or dynamically 
during the expansion or contraction of a universe or the collapse of matter to a black hole. 

One application is the behavior of universes which on larger scales has been studied from 
;hepoint of view of loop quantum cosmology by using effective classical equations 

^ly^Oil^ ■ While the fundamental description is quantum, governed by a difference equation 
for the wave function 

BBS 

, effective classical equations show the diverse cosmological 
effects more easily. The idea of using effective classical equations is that in semiclassical 
regimes they describe the position of wave packets solving the difference equation. They 
simplify the analysis considerably even in isotropic models and can be expected to do so 
even more in inhomogeneous models or the full theory. In light of the previous discussion 
an open question is where exactly an effective classical equation makes sense as a good 
approximation to the behavior of the difference equation, and where additional correction 
terms have to be taken into account. 

This question can be answered by a direct comparison of effective semiclassical descrip- 
tions, given by ordinary differential equations, with the underlying discrete quantum evo- 
lution governed by difference equations (these difference equations may even be partial 
depending on the number of matter fields). However, since an ordinary differential equation 
is quite different from a discrete difference equation, their solutions can not be compared 
directly. For such a purpose we first have to extract appropriate data from solutions of 
the difference equation, usually by taking expectation values, which we then compare to 
the classical theory or one with further corrections. (At this point one has to distinguish 
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between ambiguities resulting from choosing how to extract the semiclassical data, and out- 
right deviations from the purely classical behavior. How this can be disentangled will be 
discussed later.) In this way one can see if new effects arise or in which range one can trust 
an effective classical equation with or without certain correction terms. 

Comparing an effective theory with the underlying quantum theory is of great relevance 
to various issues in loop cosmology. For instance in order to determine the starting point 
of inflation I 111 . Il2| , observable signatures in cosmic microwave background [10| , the 

validity of effective classical bounce pictures jisl 16, III |, general aspects of the 



approach to a classical singularity |18l . Il9l | or evolution of perturbations. If the effective 
theory requires further corrections from the quantum theory which can be verified by their 
comparison, then these correction terms in principle can leave an observational signature 
which can be used to verify loop quantum cosmology. 

A particularly striking difference between the classical and the quantum theory is the 
issue of time. A common understanding which works in both cases is that of relational time, 
where time is not an external , ab solute parameter but encoded in the relative change between 



different degrees of freedom [20, |2l|, |22] . However, this concept is difficult to use explicitly, 
and so classically one employs the space-time picture where time is just a gauge coordinate. 
Thus, this time coordinate has no invariant physical meaning, but nevertheless provides a 
helpful intuitive understanding of a given gravitational system. From the Hamiltonian point 
of view, this time coordinate is the gauge parameter for orbits generated by the Hamiltonian 
constraint. In this way, coordinate time is related to the second effect of a first class classical 
constraint, namely that of generating a gauge transformation in addition to restricting fields 
to the constraint surface. 

In a canonical quantum theory the situation is different because in the Dirac procedure to 
solve first class constraints there is only one step by requiring physical states to be annihilated 
by the quantum constraints, which then are automatically gauge invariant. In systems with 
a Hamiltonian constraint, physical states are thus timeless, which has led to the name 'frozen 
formalism.' As we will discuss below, the two steps familiar from the classical procedure can 
easily be disentangled also in a quantum procedure to solve the constraints, in particular 



if the technique of group averaging 



23| is used. Physical states can then be represented in 



an evolving manner, depending on the gauge choice via the lapse function. We emphasize 
that our notion of evolution in time as used here is with respect to coordinate time and not 
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physical time. We call this coordinate time-dependent family of states a 'state-time' in order 
to indicate that classical space in a space-time has been replaced by a quantum state, while 
time remains classical. At a given time parameter of the state-time the constraints will not 
be satisfied, but the whole state-time represents a physical state, which can be reconstructed 
by integrating over time, in a well-defined way. 

As a practical application we develop a scheme to decide the domain of validity of an 
effective semiclassical description and whether in some regimes it requires further correction 
terms. For that, we compare the expectation value of, e.g., the volume in a given state-time 
with the volume obtained from the classical theory (or with an effective theory including 
further corrections). The correction terms can be derived by computing the expectation 
value of the Hamiltonian constraint operator in a coherent state and expanding around the 
classical expression 12^.125^. There are diverse sources of deviations in the case of non-linear 
constraints which can be studied analytically or numerically: First, there are the usual 
Ehrenfest statements about the relation between classical and quantum equations of motion 
for expectation values. In addition there are choices related to choosing an initial semiclas- 
sical state, and the way the constraints are violated at fixed coordinate times. Finally and 
most importantly, there are genuine quantum gravity corrections whose implications can 
give rise to new physical effects. The latter imply the modifications we want to include in 
effective classical equations. 

As we will see, it is possible to disentangle these effects at least qualitatively. The resulting 
modifications to classical behavior can then be compared to known analytical results, as done 
here for inflationary behavior and bounces, or be used to suggest new effects. In the next 
section we recall the group averaging procedure which will be followed by a brief discussion 
of the issue of introducing coordinate time in this context. In Sec. IV we will present a way 
to implement this idea numerically and study examples of universes with a cosmological 
constant or dust as matter. We will show that, both for inflation and bounces, the effective 
semiclassical theory which incorporates modified geometrical densities is in good agreement 
with discrete quantum evolution till very small scales. This proves that effects derived 
from the effective semiclassical theory capture to a large extent the true nature of quantum 
spacetime. Here we focus on describing the way our technique can prove useful in testing a 
semiclassical theory, leaving a more detailed study for future work. 
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II. GROUP AVERAGING 



We will discuss here only systems with a single constraint which we mostly think of as a 
Hamiltonian constraint generating coordinate time. In an isotropic cosmological model the 
constraint itself will be the Friedmann equation, and the evolution equation it generates is the 
Raychaudhuri equation plus equations for matter fields such as the Klein-Gordon equation 
for a scalar. After quantizing we obtain the constraint operator which annihilates physical 
states, and all evolution would have to be extracted from physical states in a relational way. 

A simple method to implement a constraint on quantum states is as follows. Consider 
the action of the constraint operator C* on a given non-physical state interpreted to promote 
a change in gauge parametrized by a real parameter A 

CM=^^M. (1) 

Obviously \(px) is not annihilated by the operator C unless \ipx) is independent of A, i.e. 
gauge invariant. To achieve A-independence we can average the state \(p) = J dX \(px) such 
that 

C\^) = j d\C\vx) = i j dA^I^A) (2) 

which in the case of compact symmetry orbits or suitable boundary conditions in the non- 
compact case vanishes identically. In conclusion, we can get a solution to our constraint by 
averaging the nonphysical state over the symmetry group: J dA Iv^a)- This is in fact not new 
because the solution to our equation for a state \ip\) is given by 

M = e-''^' \^o) (3) 

and the invariant state is 

' dXe-'^'='\ip) (4) 

which is nothing else than the group averaging map into the physical Hilbert space. 

As examples we consider different cases in which the method works with different success. 
The simplest case would be to consider the constraint Pelf) = which imposes rotation 
invariance on a two dimensional system with angular coordinate 6. Then, by following the 
prescription presented above we replace the equation 

{e\Pe\f)=tdef{e) = (5) 



with a Schrodinger hke equation. The constraint then acts on non-physical states by 

ide^x{e)=idx^x{e) (6) 

with general solution 

^x{e) = ^{\ + e). (7) 

Physical states are thus given by 

^= dA(^(A + ^)= / dA(^(A) (8) 
Jo Jo 

where any 6'-dependence is removed. 

We can also consider the non compact case by using the translation generator with 

exactly the same calculations (except that allowed functions ip have to be restricted to a 

suitable set for the A-integration to exist, similarly to selecting an appropriate subspace 

of the kinematical Hilbert space for group averaging). For a more general example let us 

consider the operator C — aXP + 6 to follow the strategy above, i.e. use the operator to 

generate transformations of an arbitrary state, and then solve the resulting partial differential 

equation. In this case the equation becomes 

(a X id^ + b) (px{x) = idxVx{x) (9) 

with solution 

^p^{x) = et(^-W-^)/(ln(x)/a + A). (10) 

Then, in order to realize the A integration we Fourier transform /(iln(a;) + A) which is 
possible only for a real: f{u) — (27r)~^ da;e~*'^"/(a;). Moreover, we commute the A- and 
a;-integrations, 

/oo /"OO J /"OO 

d\^x{x) = e'"^-^'^^^^ / "^e-'^-^'^^^^fiuj) / dAe-^(-+^)\ (11) 
■oo J -oo ^TT ^-oo 

Here, the integral in A would diverge for arbitrary complex h and we get a solution for the 
constraint only if it is real. After integrating over A and uj we obtain f^i^x) = ce*oi'^(^) which 
can be checked to satisfy our constraint. If 6 ^ R, the final integration diverges which means 
that we are not allowed to commute the integrations. Moreover, in this case we would have 
to choose appropriate fall-off conditions for /. 

Whether or not we are using a self-adjoint constraint operator has significance for the 
physical inner product, which we are not considering here. Still, adjointness properties also 
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play a role at the level of solving the constraint, as the following example given by C* = d/dx 
demonstrates. Now we have to solve the equation 

which is done by any arbitrary function ipx{x) = ip{X + ix) depending only on A + ix. 
Integrating over A to compute a physical state is now done along a line shifted vertically 
by the amount x in the complex plane, where we interpret solutions ^px{x) as holomorphic 
functions on the complex plane with coordinate z = X+ix. The result will be independent of 
X only if (f satisfies appropriate fall-off conditions and does not have residues. Both properties 
together cannot be satisfied for non-zero functions since, owing to Liouville's Theorem any 
bounded entire function is constant. Thus, any function for which the A-integrations exist 
must have poles on the complex plane such that the averaging procedure does not lead to 
constant functions of x, as expected for this constraint, but only to locally constant ones. 

While the averaging can be defined even for non- self adjoint constraints, from the numeri- 
cal point of view self-adjointness of the constraint operator is essential. Non-real eigenvalues 
would imply exponentially growing modes in solutions to the differential equation which 
lead to numerical instabilities. 



III. COORDINATE TIME 



Quantum gravity in the perspective of canonical quantization arises as a constrained 
system which gives rise to the well known Wheeler-DeWitt equation instead of a Schrodinger 
like equation. Nowadays canonical quantum gravity is usually formulated with Ashtekar 
variables and it is possible to realize the constraint algebra on a well-defined kinematical 

" — n 

1201 . The Gauss and diffeomorphism constraints can be solved 



Hilbert space 
by group averaging 



settled yet 




but discussions about the correct Hamiltonian constraint are not 
301 13 IL I32| . Also the issue of the physical inner product and how to use 
the solution space to the Hamiltonian constraint are almost completely open in the full 
theory. After reducing to symmetric models the Hamiltonian constraint simplifies and 
can often be treated explicitly. Even in the simplest cosmological models the physical inner 
product is not yet understood, but since the spatial volume can be used as internal time in 
a cosmological situation the problem of time does not play a role. In all these cases there is 
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a Hamiltonian constraint whose gauge parameter classically corresponds to coordinate time, 
and it is our aim to discuss how such a parameter can appear with this interpretation in 
quantum theory. We emphasize that, as a general problem, this is much simpler than the 
problem of time where time is understood as a physical parameter valid for the full quantum 
theory. In contrast, we are interested in formulating the quantum theory in a parameterized 
way, with a new non-physical coordinate time parameter. Moreover, this parameter is 
expected and intended to make sense only in semiclassical regimes. This will be discussed 
in more detail in the Conclusions making use of what we learned in the examples. 

As discussed before, the group averaging procedure to solve constraints can be split into 
two steps which roughly correspond to the two classical steps of restricting to the constraint 
surface and factoring out by the gauge orbits. The correspondence is not perfect, however, 
since one single member ip\^^ of a family {(y^AjAeiR of states exhibiting the gauge parameter, 
which we call state-time in the case of a gravitational system, is not a solution to the 
quantum constraint, while a classical gauge orbit can completely lie within the constraint 
surface. Even in a semiclassical regime this would lead to deviations between the classical 
and quantum behavior since the constraints are always violated when the gauge parameter 
is exhibited. Still, as an approximation and a heuristic tool the gauge dependent family of 
states can be very useful. In particular for the Hamiltonian constraint of a gravitational 
system this allows us to describe the quantum dynamics approximately (in the sense that 
the constraint is not imposed exactly for otherwise the state-time would have to be time- 
independent) with reference to a coordinate time, which justifies the use of effective classical 
equations of motion. 

Taking into account possible choices of lapse functions, which after quantizing can be 
operators if depends on t via kinematical degrees of freedom (such as N{t) = a{t) which 
is used when transforming from proper time to conformal time in an isotropic model), we 
arrive at 

rai^i) =^^|^i) (12) 

as the evolution specified by the Hamiltonian constraint operator. 

The condition that the Hamiltonian constraint has to annihilate physical states empha- 
sizes the fact that physics does not depend on (coordinate) time. Correspondingly, Eq. (fT^. 
which is directly related to a choice of time through a Schrodinger equation on non physical 
states, and the state-time |\I'i) solving it are in fact not unique because we can fix the gauge 
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freedom in different ways by making different choices of the quantized version N of the lapse 
function. 

If we were interested in completing the group averaging we would have to integrate the 
state-time over t in order to arrive at a physical state. However, just as the classical space- 
time picture, which we are interested in here, arises only when gauge orbits parameterized 
by coordinate time t are not factored out, we have to refrain from performing this final 
step and instead work only with the state-time. In principle, with much effort one can 
always remove all gauge dependence, classically by factoring out gauge and in quantum 
theory by integrating over the gauge group, but physical intuition is best developed in a 
(coordinate) time dependent picture. Nevertheless, we think that the discussion of group 
averaging justifies our use of the state-time and the following applications. 



IV. APPLICATIONS 



For semiclassical physics it is very convenient to have an explicit coordinate time pa- 
rameter since classical intuition is based on the space-time picture. In principle it is also 
possible to work with internal times both at the quantum and classical level, but it comes 
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in fact, most 
are based 



with much more technical effort. In the case of loo 
recent phenomenological applications 
on effective classical equations which are differential equations in coordinate time 

and implement the main non-perturbative quantum effect in matter Hamiltonians 



Direct studies of the underlying difference equations, on the other hand, are more 



complicated 
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These equations show the main qualitative effects that have to be expected from the loop 
quantization at a more intuitive level. For more quantitative applications, however, it is 
important to see to which degree they provide an approximation (in the sense discussed in the 
Introduction) to the quantum behavior and whether in some regimes additional correction 
terms have to be taken into account. One way of evaluating these correction terms for 
instance is to compute the expectation value of the constraint operator in a coherent state 



25|. This results in the classical constraint together with correction terms, and thus corrected 



equations of motion for, e.g., the scale factor. The present paper suggests an alternative 
procedure, which is one step closer to the quantum theory. Using the quantum coordinate 
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time picture, we can evolve the state first within the quantum theory, and then compute 
the expectation value of, e.g., the volume operator which also gives us the time dependence 
of the scale factor. Both procedures are approximations to the quantum dynamics: In the 
first case one uses kinematical coherent states to compute the expectation value, while in the 
second case coordinate time is introduced which, as explained above, is not exact in quantum 
theory. In both cases, however, the Hamiltonian constraint is used and implemented at least 
partially: The first case imposes a classical constraint with quantum corrections, while in 
the second the quantum constraint is used to evolve states, which then could be averaged if 
we are interested in the physical state. 

In a sense, the procedures differ by a commutation of evolving and translating from quan- 
tum to classical behavior (by taking expectation values). We either take the expectation 
value of the constraint first (in a kinematical coherent state) and then determine the evolu- 
tion in a classical manner, or we first evolve a kinematical state with quantum operators and 
then determine classical quantities from expectation values in the resulting, time-dependent 
states. It is not guaranteed that the two different steps in fact commute, which leads to dif- 
ferences between the two procedures. In simple models, however, both ways of determining 
the dynamics can be implemented at least numerically and then compared with each other. 
We will demonstrate this in what follows, leaving a more detailed investigation with precise 
statements of ranges of applicability and of the necessity of additional correction terms for 
future work. 



A. Numerical implementation 

In isotropic loop quantum cosmology for a flat model, the Hamiltonian constraint operator 
is given by 0, Q 



(^oV^)m = - v^+3)v^M+4 - 2(v;.+i - v^-i)v^^ 



+ (V"^_3 - V^-5)^A'-4 + |7rG'7^4 Sgn(/i)i^matter(/i)^A' 



acting on a wave function ■?/': M — > C supported on eigenspace of the triad operator p. 
The coefficients are given in terms of the volume eigenvalues = (i7£p)^/^|/ip/^, with the 
Barbero-Immirzi parameter 7 = 0.238 



entropy 



44 



42 



as it follows from calculations of black hole 



and iymatter(/^) IS the matter Hamiltonian which we will choose later. 
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Since only labels with distances of four apart are involved, we introduce 4m := /i with 
integer m (see also 0] for this restriction) and write the operator as 

{Holl!)m = {Vm+5/4 — Vm+3/4)'^m+l — '2{Vm+l/4 — Vm-l/4:)'^m 

+ {Vm-3/4: - Vm-5/i)i^m-l + |vrG7^4 Sgn(m)i^matter (4m)?/'m 

with ipm '■= V'4m ^iid Vm '■= V4^rn- This Operator will have to be symmetrized to H := 
^{Hq + Hq) for the numerical implementation of coordinate time to be stable. (Alternatively, 
we can choose a lapse function N{t) = a{t), corresponding to conformal time, and quantize 
such that N has eigenvalues proportional to V^+i — V^-i. The resulting operator NHq would 
be symmetric without reordering Hq^ but since the lapse function and its quantization vanish 
for a = 0, the quantum behavior around /i = will be problematic.) 

For numerical purposes we need to restrict the operator to finite lattices C^^^m '■= {n & 
Z : mc — N/2 < m < mc + N/2} of size and centered at rric, such that it will be represented 
by a tridiagonal x A^-matrix Hij with Hij = for j > i + 1 or z > j + 1, 

Ha = -2{Vm,-N/2+i+l/4: - Vm,-N/2+i-l/4) + |vrG7^4 Sgn(mc -^N + i)-ffmatter("^c - + i) 

(13) 

and 

= Hi^i^i = \{Vmc-N/2+i+5/4. ~ ^mc-iV/2+i+3/4 + Wmc-Af/2+j+l/4 ~ ^^mc-Af/2+i-l/4) • (14) 

The above splitting of Hamiltonain matrix into diagonal and off-diagonal parts is done purely 
for numerical convenience. We now introduce the coordinate time parameter t and solve 

H-^t = i^^tPt (15) 

numerically. The solution to the differential equation (fT3j) is then given by 

^pt = exp{-itH) ■ V'in 

with an initial state tpin G C^. The function tpt can then be used to compute time dependent 
expectation values. 

B. Examples 

As examples of our technique, we would now investigate various cases of matter Hamil- 
tonian and compare the effective semiclassical theory with the evolution determined by 
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quantum difference equations. We first consider the case where matter Hamiltonian is just 
a cosmological constant and the scale factor is the sole degree of freedom. Unlike with an 
internal time evolution, for which no other degree of freedom besides internal time a would 
be left, this still allows us to have a non-trivial evolution of a in coordinate time t. This 
example is then followed by examples of inflation and bounce with inclusion of matter (dust) 
in the analysis which signifies the role played by effective densities in a good semiclassical 
description. 

The total Hamiltonian is given by 

^ad2-87rG'i/^atter = (16) 



which leads to the Friedmann equation 



p{a) . (17) 



0? 3 

In our first example we take p{a) = A and later we will consider the classical dust density 
p{a) — Mjc? . In general, there are two qualitatively different kinds of modifications which 
bring the behavior of effective cosmological equations closer to that of the quantum theory, 
which we will study in what follows. First, by non-perturbative effects the geometrical 
density in the matter Hamiltonian density is replaced by a function which is finite and 
docs not diverge at small a. The precise form of the function will be discussed below. 
Secondly, there are perturbative corrections in the gravitational part d^/a^ of the constraint 
which appear as additional terms on the left hand side of the Friedmann equation. We will 
mainly discuss one such term, which includes effects of the spread of a wave packet in the 
effective classical framework. 



1. Effective classical behavior: Cosmological constant 

In order to compare the evolution given by an effective theory and the quantum difference 
equation we first have to choose an initial state ipi^. This state should be peaked on a specified 
classical volume, together with an extrinsic curvature which follows from the volume and the 
classical Hamiltonian constraint to which we compare the quantum evolution. In principle, 
one can use any suitable state as an initial state, however for simplicity we will use a Gaussian 

V'in.fe = exp ^-i^^^^ + 2i{k - N/2)co 
12 



in what follows. Here the isotropic connection component cq is related to extrinsic curvature 
by c ^ ^.an. It .s determined fro,„ the peak seale faetor a„ ^ ^/|^ M by the 

constraint —Qc^ao + 87rG7^ifmatter(ao) = 0. (At this point the classical equations to be 
compared with enter. If correction terms are included, as will be done later, the value of 
Co changes and so does the quantum evolution of the new initial state.) Nevertheless, the 
explicit choice has an influence on the evolution such as the degree of spreading of the wave 
packet. Numerically, a choice leading to less spreading will allow the evolution to be reliable 
for a longer period of time since the boundary values will remain small longer and finite size 
effects will set in later. 

For our purpose, comparing the quantum evolution with classical equations, the choice 
of initial wave packet is not that important since it is sufficient to know the evolution for 
limited amounts of time only. We will not be able to compare whole solutions as functions 
of time in this way, but since the classical equations and also the corrected ones are local in 
time we can study deviations. One can then decide which correction terms are needed to 
describe the local change in volume following from the quantum evolution. 

A quantitative statement about the deviations is complicated by the fact that the quan- 
tum evolution is not uniquely related to a classical expression. One can take expectation 
values and compare with the classical functions, but due to the spread of the probability 
distribution the result depends on whether we take, e.g., the expectation value 

N 
k=l 

of the volume operator, or that of the scale factor operator a = V^^^ and the cubic power 
afterwards. 

Using these ingredients we compute the volume for the quantum equation (|TH|l and the 
classical theory feqll7|). Classically, we have V = a^, but this relation will certainly not be 
satisfied for the expectation values of V and d once the spread of the wave function becomes 
large. In Fig. ^ we have shown the behavior of expectation values of the volume operator 
V (+) and that evaluated from d (x). These are compared with the volume calculated from 
the classical theory (dashed curve). As it is clear, the classical theory gives a very good 
approximation to the underlying evolution from the difference equation. However, due to 
spread of the wavepacket some discrepancy occurs for late times. This can also be noticed 
in the evolution of the wave packet which is shown in Fig. |21 Comparing the two different 
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FIG. 1: Expectation value of the volume (+) compared to the effective solution V{t) (dashed) 
and the volume (a)^(i) computed from the expectation value of the scale factor (x) for the case 
of cosmological constant. The initial peak of the coherent state is chosen around rric = 200 with 
(7 = 20 and iV = 500. We take A = 10-^. 



ways of computing volume expectation values shows that the discrepancies do not result 
from new physical effects but only reflect the ambiguous way of relating expectation values 
to classical behavior. 

In general, departures between expectation values of the volume operator and classical 
values can have several reasons. Besides approximations used in the method to solve the 
constraint, there are quantum effects which we are most interested in here. They can be 
divided into two classes, the first one arising from small-scale or high-curvature effects in 
quantum operators, the second one coming from the fact that we have an evolving wave 
packet rather than a sharp classical point. Both effects can be included into effective classical 
equations, but in the second case it is not always clear if modifications to the classical 
behavior are a consequence of having chosen a bad initial state or a physical effect related 
to the evolution of spread, asymmetry and other deformations in profile of the probability 
distribution. 
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FIG. 2: A spreading wave packet for A = 10~^ as in Fig. ^ At the back of the right hand side 
one can see small oscillations building up when the wave reaches the boundary. The vertical axis 
shows the magnitude of the normalized wave function, which starts as a Gaussian at the front and 
evolves to the back. 

In order to quantify this, one can evaluate the skewness of the wavefunction as it evolves. 
The skewness which is initially zero for the Gaussian coherent state describes the asymmetry 
of the wave packet and is given in terms of various expectation values of powers oi p = V"^^^ 
as 



a3 



(p')-3(p)(p2)+2(p)3 (19) 

where a is the standard deviation ((p^) — (p)^)^^^ which initially is given by a^^ = |7cr. (We 
use p in order to define skewness because it is the isotropic component of the densitized 
triad, which is basic in loop quantum gravity. Eigenvalues of p are proportional to /i or m 
such that the skewness is computed for the variable in the wave function.) It is clear from 
Fig. El that the initial Gaussian wave packet gets skewed with time and deforms, but the 
deformation remains small such that the main parameters characterizing the wave packet in 
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FIG. 3: Growth of standard deviation {a /din) and skewness (s) for the cosmological constant case 
studied in Fig. ^ The coherent state remains almost peaked over the classical value due to small 
growth in skewness. 

the cases studied here are only its expectation value and spread. 

The evolution of skewness with time reflects the fact that the coherent state is no longer 
peaked at its classical value, which is one cause of departure between quantum and classical 
curves. Fig. |21 shows the spread of the wavefunction with respect to its initial value. A 
symmetric spreading would give an evolved coherent state still peaking at classical values 
at late times for the operator a? corresponding to the discrete argument m of the wave 
function. The signiflcance of small growth in skewness for the time scale of interest is that 
the probability distribution of the wavefunction remains peaked close to the classical volume. 
Moreover, since the quantized volume and scale factor are given by powers of m different 
from one, they do not follow the classical curve exactly, their difference increasing with 
increasing spread and skewness. 

This discussion shows that a comparison between quantum and (effective) classical be- 
havior cannot be done arbitrarily precisely because of the unsharp nature of quantum wave 
packets. To determine the level up to which a comparison is reliable we have different tech- 
niques as illustrated in this case. The spread of the wave packet and its deformation can 
be computed and plotted such that a strong growth signals stronger departures. Similarly, 
plotting expectation values of different powers shows a window in which the classical curve 
has to be expected if there would be no quantum modiflcations to the equations of motion. 
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As we will see in the following examples, effects from such quantum modifications are much 
stronger such that they can easily be separated from simple spread effects. 



2. Accelerated expansion 



In the previous subsection on cosmological constant we showed that evolution extracted 
from the quantum difference equation by introduction of coordinate time agrees well with 
the evolution determined classically. Now we will study the case of dust Hamiltonian by 
introducing the modification to geometrical density in the Friedmann equation ()17p . This 
modification to geometrical density is a novel prediction of loop quantum cosmology at shor ; 



scales and has led to various interesting applications 

llfil . it! lis! loj l . The main new ingredient is then given by a non-perturbative modification 
to the classically diverging density (which has been derived in 3^ and is based on 
expressions of the full theory 3^) relevant at small volumes. If a matter Hamiltonian has 
a density term then such modifications would enter dynamics and change the behavior at 
small scale factors. For example, the Hamiltonian for a massive scalar field is given by 

H^ia) = ^a-'pl + a'Vi<l)) (20) 

with momentum and potential V{(f)). This Hamiltonian in loop quantum cosmology is 
modified to 



where dj^ is the modified geometrical density given by 



(21) 



(22) 



where 



A(g) 



2/ V/ + 2 



[(9 



\l+2 



sgn(g-l)|g-l|'+i] 



3/(2-2/) 



(23) 



and j is a quantization ambiguity parameter (a half- integer) . There is another ambiguity 
parameter < / < 1 [2^ , which is more restricted by full loop quantum gravity and is usually 
taken as / = 3/4. For very small a -C y/jip, dj^i{a) behaves as a positive power of a, 

3 



dj^i{a) ~ 



3/(2-20 



v3(2-0/(2-20^-3 



(24) 
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I.e. 

d,{a) := d^. 3(a) ~ (y)' (|74j)-^'/'«^' • (25) 

This is the main ingredient for effective classical equations with the most dramatic effects, 
and the methods developed here allow us to put this term in effective equations on a more 
solid footing. In order to avoid numerical complications related to the additional degree of 
freedom (p we use a dust model and employ the same effective density for its energy density 
p{a) = M/a^. We emphasize that for dust this modification to the matter density is not the 
one most naturally expected from loop quantum gravity. Rather, in loop quantum gravity 
the matter Hamiltonian is primary and will be quantized. Since this is a constant M for 
dust, there would be no modifications of this kind at all to the classical equations. We use 
the modification (which can also be interpreted as arising from an additional quantization 



ambiguity 







47l |) here to model the kinetic term of a scalar field Hamiltonian, and in order 
to study the implications of effective densities. 

To that end we replace the classical matter energy density for dust, p{a) = M/a? with 
a constant M, by p(a) = Mdj i{a) and study the evolution as determined by (fT3j) and the 
effective theory. Note that unlike the case with a cosmo logical constant, the effective theory 
is different from the classical one (in which p{a) = M/a^) due to the introduction of a 
modified geometrical density. To see the differences between the two (which also highlights 
the role of introducing dj^) let us first compare the volume expectation values from ()15p with 
the classical theory. We have shown the results in Fig. ^ It is interesting to see the way 
the effective density modifies the quantum dynamics given by difference equations. As is 
clear from plot, the classical theory does not match the quantum description which indicates 
that the volume expectation value increases more strongly than the classical solution, which 
unlike in Fig. ^ is not an effect of a spreading wave packet or skewness. Fig. shows that 
spreading and skewness of the wave packet are not large. This is also depicted in Fig. IHl 
showing that wave packet does not spread strongly during the displayed evolution. 

From Fig. 13 it can be seen that the expectation value of the volume operator agrees better 
with the cubic power of the expectation value of the scale factor than with the classical 
solution. This is another sign that the modified dynamics is responsible for the deviations, 
rather than just change in shape of wave packet during evolution. However, Fig. |3] does 
not tell us decisively how well the expectation values would agree with modified classical 
dynamics. A hallmark of the effective density dj{a) is that it implies inflation (accelerated 
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FIG. 4: Expectation value of the volume (+) compared to the classical solution V{t) (dashed) and 
the volume computed from the expectation value of the triad (x) for dust with M = 10 and an 
initial peak around m-c = 100 (A^ = 1000, a = 20). The ambiguity parameter j for the effective 
density is j = 400, such that the density peaks at w j/2 = 200 corresponding to a volume 
K = (27^pm*/3)'^/^ « ISOip, which is reached around t = 1.1. 

expansion) when inserted into the classical equations and when the scale factor is below the 
peak value. In fact, the expectation values in Fig. E] increase more strongly than the classical 
solution, but in the range < t < 1.1, where we are below the peak, deviations are small. 

It is much more illuminating to plot time derivatives, or rather difference quotients from 
the numerically obtained data, and compare with the classical behavior. Fig. [7| shows that 
in fact the derivative of the scale factor increases when t < 1.1, i.e. below the peak of dj{a), 
and thus accelerates in agreement with the expectation. This figure also shows that classical 
description completely fails to capture the variation of scale factor with time as dictated by 
quantum theory. 

So far we have seen that the quantum evolution differs significantly from the classical 
one and in particular leads to accelerating scale factor expectation values and thus inflation. 
Now we turn our attention to the effective theory where in the matter density is replaced 
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FIG. 5: The evolution of standard deviation and skewness for the case in Fig. ^ The spread 
increases most strongly in the inflationary regime (t < 1.1)), while the wave packet is skewed 
stronger during unaccelerated expansion. 




FIG. 6: A spreading wave packet for the evolution shown in Fig. [l] 
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FIG. 7: Time derivative of scale factors (+ corresponds to those computed from volume expectation 
value and x refers to those of scale factor expectation values) plotted in Fig.0J Dashed curve shows 
the classical curve. The classical description does not match the evolution from quantum theory. 

by dj{a). We thus obtain a modified classical description which as we will show agrees quite 
well with evolution determined by difference equations. 

In order to check the effective density more directly we compare the expectation values 
to numerical solutions of this effective classical equation with the matter density replaced 
by Mdj{a). The result is shown in Fig. |H1 On comparison with Fig. IH it is clear that the 
effective theory matches the quantum evolution much better than the classical description. 
We have plotted the derivative of scale factor with time in Fig. El which shows that time 
variation of scale factor as computed from effective theory agrees well with the change of the 
expectation values of scale factor in the regimes before and after the peak. In particular, the 
inflationary behavior in the modified region with increasing a can be seen from the effective 
semiclassical as well as the quantum solution. After the peak, both show the expected 
non- inflationary behavior. 

Interestingly, around the peak the effective classical time derivative is larger than the 
change in expectation values, which is also the reason why the effective classical volume is 
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FIG. 8: Volume expectation values (+ for {V) and x for {a)'^) compared to the effective classical 
solution (dashed). Compared to the points in Fig. ^ the quantum behavior is different since the 
initial wave packet is now peaked on the effective rather than unmodified classical constraint surface 
(which decreases the initial cq). 

slightly larger than the quantum volume at later times in Fig. |H1 At the point where the peak 
occurs, d is largest so that higher order corrections (higher powers of a in the Friedmann 
equation) are expected to have the strongest influence. Those corrections have not yet been 
included into effective semiclassical equations and we leave a more detailed investigation of 
their effect around the peak for future work. Nevertheless, one can see numerically that in 
the case studied here the effect of higher order terms is negligible. This can also be seen 
from the fact that those higher order terms arise as powers in c = —^70 which thanks 
to the smallness of 7 remains sufficiently small compared to one throughout the evolution. 
Moreover, for the flat model the Hamiltonian is invariant under change of sign in c such that 
the next higher order correction is suppressed by a power c^. 

Rather, the reduced derivative of the scale factor around the peak is a consequence of 
deformations of the wave packet as shown in Fig. ^1 In fact, around the peak the effective 
density changes rather rapidly from increasing to decreasing behavior. Thus, a part of the 
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FIG. 9: Time derivative of scale factors plotted in Fig. |H1 Inset shows variation of Hubble rate 
with time {H = a /a). 




FIG. 10: The time variation of standard deviation and skewness for dust with effective density 
modification. 
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wave packet will already be in the decreasing regime while its center is still in the increasing 
one, which lowers the overall density seen by the wave. This behavior is verified by looking 
at the skewness of the wave packet during the evolution which is plotted in Fig.^1 First, the 
skewness turns positive which means that the right tail of the wave packet becomes heavier 
than the left one. At some point, skewness starts to decrease and reaches negative values, 
describing a redistribution of parts of the wave packet such that now the left tail becomes 
more pronounced. This redistribution implies that expectation values of powers of m, such 
as the scale factor and volume, are lowered as compared to the expected evolution. In fact, 
comparing Fig. El with Fig. IHl shows that the turn-around of skewness starts just when the 
deviations between the change of (a) and the effective a appear, with the expectation values 
increasing less strongly than the effective classical value. 

This observation shows that effective classical equations used so far do not capture all 
details in the quantum evolution around the peak. On the other hand, the behavior before 
and after the peak is described very well. Cosmological studies so far have mainly focused 
on the modified behavior at small scale factors and the initial inflationary epoch, which 
is modeled reliably by using just the effective density. The peak behavior was actually 
more problematic since the Hubble parameter in the case of scalar dynamics easily became 
dangerously high, larger than one in terms of Planck units, which leads to doubt in the 
further semi-classical evolution. Here, the quantum behavior with smaller a suggests that 
additional effects from the wave packet can lead to a better semiclassical picture, which may 
be modeled in effective classical equations by taking into account effects of skewing wave 
packets. How this appears in the case of a scalar field, and whether in this case a can become 
large enough for higher order corrections to be relevant, remains to be studied. 

3. Bounces 

So far we only looked at the evolution for rather large volume. For smaller |m|, the 
approximation by classical evolution will become worse and eventually break down. When 
exactly this is happening depends on the parameters for the cosmological model and the 
choice of initial wave function. As an example we again use the dust model but evolve to 
earlier times, towards the classical singularity. In Fig. ^2 we show the evolution of the wave 
packet as it deforms once a significant part of it reaches m = 0. As in the case of dust we 



24 



FIG. 11: Wave packet evolving toward the classical singularity and bouncing off, only penetrating 
negligibly to negative m (to the right). The parameters are rric = 100, N = 500, o" = 10, M = 10, 
j = 200. 



first compare volume expectation values obtained from ()15p with the classical theory. The 
result is shown in Fig. which shows a bounce at non zero volume for expectation values. 
The classical curve first hits the singularity at zero volume. 

It should be noted that the effective classical equations of a fiat model do not allow a 
bounce even when we use the effective density in the dust case, and thus this bounce is not of 

gnnnPi 



the semiclassical type as in 



14 Il5l lid [iTf . Further, deformations of the wave function 



show directly that the classical space-time picture is not valid during the bounce phase. 
Moreover, the nature of a wave packet as opposed to sharp classical values becomes more 
relevant, which provides a quantitative explanation of the observed bounce. The spread a 
implies a correction to the Friedmann equation which then takes the form 



. 2 I 1 -2 -2 

a + ^7 a 



(26) 



where for simplicity we ignore the a-dependence of a due to spreading. Using the small-a 
expansion for dj{a) in p(a) = Mdj{a) and setting d = for the bounce results in a 
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FIG. 12: Volume expectation values for the wave packet in Fig. 111! The expectation values bounce 
at non-zero values, while the classical curve hits zero volume. 

bouncing scale factor (in Planck units) 

abounce ~ 3-2^/2« {rf' ila-'M-V'^Y'^'Y^'' ^ 2.44. (27) 

Since we used the small-a expansion, which overestimates dj{a), the bounce value is slightly 
larger resulting in Obounce ~ 2.50 if the function dj{a) is used in its full form. The correspond- 
ing bounce volume, Vbounco ~ 15.6 is considerably smaller than the minimum expectation 
value in Fig. ^| But if we use the effective Friedmann equation fl26|) to place the initial 
wave packet on the effective constraint surface the expected bounce radius and the numer- 
ical one in Fig. IT!^ agree within the limits of the expectation values (V^) and (a)^. In fact, 
the wave packet bounces earlier than expected, but around the expected volume. Thus, the 
spread-dependent correction term explains why there is a bounce and gives a good estimate 
for the bounce radius. 

Nevertheless, it is also clear from Fig. ^1 that the agreement between the expectation 
values and the effective solution deteriorates around the bounce. This can also be seen from 
Fig. El which shows considerable change in spread of the wavefunction. In fact, the wave 

26 



X i ± * 




FIG. 13: Volume expectation values for the wave packet initially peaked at the effective constraint 
surface compared to a numerical solution (solid line) of the effective Friedmann equation (|26() . 
Since the corresponding wave packet, plotted in Fig. 1151 moves closer to the classical singularity 
than in Fig. IllL the spread of wavefunction increases faster than in Fig. 1121 
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FIG. 14: Increase in spread and skewness for backward evolution through a bounce as in Fig. 1131 
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FIG. 15: Wave packet evolving toward the classical singularity and bouncing off as in Fig. 1111 but 
initially peaked on the effective constraint surface at nic = 80. Since the bounce radius now is 
smaller, the leakage to negative m is larger (see rightmost contour line). 

packet does not only spread but also separates into different packets as seen in Fig. UHl Thus, 
even taking into account the spread dependence in the modified Friedmann equation would 
not completely describe the quantum behavior. (Note that the spread dependent correction 
term was derived under the assumption of a Gaussian wave packet. Fig. also shows that 
the skewness does not increase strongly, but for this case of a wave separating into different 
packets skewness alone is not a good measure for the deviations from a Gaussian.) The 
evolution can then no longer be seen as semiclassical and it is not possible to get a better 
agreement by including further corrections into an effective equation. 

Still, one can also see that there is a rather undisturbed wave packet bouncing off, while 
other parts of the wave function stay around the classical singularity which does not affect 
expectation values of geometrical quantities very much. Thus during the evolution shown 
here the expectation value of the volume and the cube of that of the scale factor do not 
deviate too much near the classical singularity. However, since strong oscillations build 
up rapidly, there are strong curvature fluctuations. This is shown in Fig. ^1 where the 
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FIG. 16: Extrinsic curvature fluctuation Aa = y ((a)^) — (ct)^ corresponding to Fig. 1131 where the 
operator for a is derived from the c-operator mapping a state ipm to |i(V'm+i — V'm-i)- 



fluctuations of a, computed from the operator 



(28) 



which initially are smaller than a increase to have a maximum at the bounce. After the 
bounce, curvature fluctuations decrease but stay larger than initially, and are comparable 
to the value of a in Fig. [T7| We thus can still think of semiclassical spatial slices of a certain 
volume at least in early stages of the bounce, but since the extrinsic curvature is not sharp 
one cannot think of them as forming a classical, smooth space-time. 

If, as before, we compute the time derivative of the expectation value of the scale factor, 
rather than the expectation value of curvature, it is seen from Fig.Elthat it is still increasing. 
That is it accelerates, as a consequence of the effective density. This figure also shows that 
the expectation value of the d-operator ()28|) follows the change of (d) more closely than 
the effective classical solution, which implies that the interpretation of extrinsic curvature 
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FIG. 17: Time derivative of the scale factor expectation values in Fig. compared to the effective 
solution (solid). The dashed line connects expectation values (a). 

(computed from the geometrical quantity a as compared to connection components) has an 
unambiguous meaning. 

As seen in Fig. the wave function leaks only slightly into the part of minisuperspace 
corresponding to negative m (i.e. the part corresponding to the side beyond the classical 
singularity in the internal time formulation ^ l^^). This is another reason for the fact that 
we can still approximately speak of a classical bounce of the volume. If the wave function 
at negative m would be larger, also the spread in spatial geometrical quantities and not just 
in extrinsic curvature would be large and classical geometry would break down completely. 
Whether or not this is happening depends very sensitively on detailed quantization choices 
in the Hamiltonian constraint and initial values of the wave function (see also 38, 41, 3|). 



V. CONCLUSIONS 



Splitting group averaging into two steps, first evolving an initial state and then integrating 
along the group orbit, allows us to introduce a coordinate time parameter into quantum 
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gravity, although only in an approximate sense. The resulting time-dependent family of 
states, the state-time, does not solve the constraint exactly; only integrating over time 
results in a physical state. Moreover, this does not introduce a physical time parameter into 
the quantum theory, but only coordinate time in semiclassical regime. Still, the resulting 
evolution equations are helpful in semiclassical analysis and in providing intuitive pictures 
of quantum effects. 

Our main application in this paper is a justification of using effective densities as the 
most prominent quantum gravity effect imported into effective classical equations. As it 
turns out, the modified classical equations describe the quantum evolution of wave packets 
very well, down to surprisingly small scales. Stronger deviations become noticeable only 
when the wave packet starts to touch the classical singularity. When exactly this happens 
depends on details of the model and the chosen initial state. The diverse effects giving rise 
to departures from classical behavior can be separated by varying the parameters involved 
in the models. For instance, we chose rather large values for the ambiguity parameter j 
determining the peak of the effective density, in order to distinguish this effect from that 
of perturbative corrections. In this way it is possible to study each correction term in the 
effective Hamiltonian separately. 

The replacement of classically diverging factors of in dynamical equations by a 
bounded effective density dj^i{a) has been confirmed as the main effect. Effective densi- 
ties themselves are partly responsible for this observation since, via the effective Priedmann 
equation, they lead to a reduction of extrinsic curvature a. This means that higher order 
corrections, i.e. higher powers of a, are less relevant at small a than without effective den- 
sities. Visible effects then occur most likely around the peak of the effective density, where 
a is largest. As observed here, the quantum Hubble rate is in fact smaller than would be 
expected just from the effective density, which is helpful for cosmological applications. How- 
ever, numerically one can check that in the cases studied here the reduced Hubble rate is a 
consequence of deformations of the wave packet, parameterized here by the skewness, rather 
than of higher order terms. Such an effect would have to be included in effective classical 
equations by introducing a correction term depending on the skewness of an evolving wave 
packet. That it is possible to describe the influence of properties of wave packets on the 
evolution by effective classical equations has been demonstrated here by studying bounces 
implied by a spread dependent correction term. The corresponding correction for skewness. 
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however, is not known since in the derivations of correction terms so far Gaussian states 
with zero skewness have been used j^]. For cosmological purposes it would be interesting 
to apply this technique to the case of a scalar field and study the role played by effective 
densities in particular around the peak. 

Detailed formulas for perturbative correction terms, which include higher order and higher 
derivative corrections, and uncertainty correction terms related to the spread of the wave 
function, are currently being evaluated. Once available, they can be used for a direct compar- 
ison with the quantum evolution as studied here. In particular the uncertainty corrections, 
which are proportional to cr^^ [25], play a role at small volume since, unlike higher order 
corrections, they are not suppressed by effective densities. An application of this correction 
term can be found in studying bounces, where it provides an explanation for the bounce 
in Fig. ^2 which does not follow from effective densities alone. The expression of the 
bounce radius determines the scale where quantum effects from the wave packet become 
important, indicated by correction terms depending on the spread a. For a phenomenolog- 
ical analysis, flbounce can be used as initial value of the scale factor, which is relevant for 
estimates of the amount of inflation. Since the value derived here for dust depends on all 
the parameters of the model, in particular the ambiguity parameter j, it is clear that the 
analysis will be more complicated than the original ones, where the initial scale factor was 
assumed to be ciinitiai = \fl^v- With an expression like Obounce, niore reliable estimates can 
be made and constraints on the parameters can be tightened. Thanks to the strong initial 
increase of dia) with a, which is responsible for the small power 1/14 in (j2Z|), the dependence 
on parameters of the model is, fortunately, rather weak. 

At very small scales the evolution of wave packets shows when the classical space-time 
picture can be trusted, when it needs to be corrected and when it breaks down completely. 
In this regime the evolution becomes much more sensitive to quantization choices in the 
Hamiltonian constraint. In particular this is true for the leakage into the domain of negative 
m, where the orientation of space is reversed, which is not surprising since this corresponds 
to evolution beyond the classical singularity in the internal time picture. In this context one 
should note that we had to use a symmetric ordering of the constraint for the coordinate 
time evolution to be numerically stable. This ordering already implies changes to the issue 

and on the relation between the wave function at positive and at 



of initial conditions 
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negative m. The alternative procedure of using a lapse function Nit) = a{t) and quantizing 
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NH symmetrically, as mentioned in Sec. IIV Al also changes the issue of initial conditions 
since the lapse vanishes at the classical singularity and the wave function at negative m 
completely decouples from that at positive m. Thus, while the coordinate time picture is 
well suited to justifying effective classical equations at non-vanishing volume, the issue of 
the classical singularity can be understood only by using the wave function directly and thus 
employing an internal time to formulate evolution. Since the classical space-time picture 
breaks down in this regime, there is no analog to coordinate time. 
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